home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dlasq2.z / dlasq2
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAASSSSQQQQ2222((((3333FFFF))))                                                          DDDDLLLLAAAASSSSQQQQ2222((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLASQ2 - DLASQ2 computes the singular values of a real N-by-N unreduced
  10.      bidiagonal matrix with squared diagonal elements in Q and  squared off-
  11.      diagonal elements in E
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DLASQ2( M, Q, E, QQ, EE, EPS, TOL2, SMALL2, SUP, KEND, INFO )
  15.  
  16.          INTEGER        INFO, KEND, M
  17.  
  18.          DOUBLE         PRECISION EPS, SMALL2, SUP, TOL2
  19.  
  20.          DOUBLE         PRECISION E( * ), EE( * ), Q( * ), QQ( * )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.         DLASQ2 computes the singular values of a real N-by-N unreduced
  24.         bidiagonal matrix with squared diagonal elements in Q and
  25.         squared off-diagonal elements in E. The singular values are
  26.         computed to relative accuracy TOL, barring over/underflow or
  27.         denormalization.
  28.  
  29.  
  30. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  31.      M       (input) INTEGER
  32.              The number of rows and columns in the matrix. M >= 0.
  33.  
  34.      Q       (output) DOUBLE PRECISION array, dimension (M)
  35.              On normal exit, contains the squared singular values.
  36.  
  37.      E       (workspace) DOUBLE PRECISION array, dimension (M)
  38.  
  39.      QQ      (input/output) DOUBLE PRECISION array, dimension (M)
  40.              On entry, QQ contains the squared diagonal elements of the
  41.              bidiagonal matrix whose SVD is desired.  On exit, QQ is
  42.              overwritten.
  43.  
  44.      EE      (input/output) DOUBLE PRECISION array, dimension (M)
  45.              On entry, EE(1:N-1) contains the squared off-diagonal elements of
  46.              the bidiagonal matrix whose SVD is desired.  On exit, EE is
  47.              overwritten.
  48.  
  49.      EPS     (input) DOUBLE PRECISION
  50.              Machine epsilon.
  51.  
  52.      TOL2    (input) DOUBLE PRECISION
  53.              Desired relative accuracy of computed eigenvalues as defined in
  54.              DLASQ1.
  55.  
  56.      SMALL2  (input) DOUBLE PRECISION
  57.              A threshold value as defined in DLASQ1.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAASSSSQQQQ2222((((3333FFFF))))                                                          DDDDLLLLAAAASSSSQQQQ2222((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      SUP     (input/output) DOUBLE PRECISION
  75.              Upper bound for the smallest eigenvalue.
  76.  
  77.      KEND    (input/output) INTEGER
  78.              Index where minimum d occurs.
  79.  
  80.      INFO    (output) INTEGER
  81.              = 0:  successful exit
  82.              < 0:  if INFO = -i, the i-th argument had an illegal value
  83.              > 0:  if INFO = i, the algorithm did not converge;  i specifies
  84.              how many superdiagonals did not converge.
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.